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Abstract. We introduce a Sherrington-Kirkpatrick spin-glass model with the 
addition of elastic degrees of freedom. The problem is formulated in terms of an 
effective four-spin Hamiltonian in the pressure ensemble, which can be treated by 
the replica method. In the replica-symmetric approximation, we analyze the pressure- 
temperature phase diagram, and obtain expressions for the critical boundaries between 
the disordered and the ordered (spin-glass and ferromagnetic) phases. The second- 
order para- ferromagnetic border ends at a tricritical point, beyond which the transition 
becomes discontinuous. We use these results to make contact with the temperature- 
concentration phase diagrams of mixtures of hydrogen- bonded crystals. 
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1. Introduction 

Orientational models with disorder and elastic degrees of freedom have been used to 
investigate phase transitions in a number of systems, as the mixed molecular crystal 
K(CN)^Bri_a; [ij. The work on quadrupolar glass models in the presence of random 
strain fields, which are supposed to mimic the random mismatching of molecular 
groups, was reviewed by Binder and Reger [2]. A compressible spherical model, with 
random bonds and in random fields, has been introduced to account for the peculiar 
behavior of compositionally disordered perovskites, also known as relaxor ferroelectrics 
[3]. Disordered Ising models, with random competing interactions and elastic degrees of 
freedom, have also been used to account for the phase diagram and the glassy transition 
in mixtures of ferro and antiferroelectric hydrogen-bonded crystals of the KDP family 
[HEIE]. More recently, random quadrupolar models in the presence of anisotropic strain 
fields have been studied in the context of nematic liquid crystal elastomers [7]. These 
investigations on disordered compressible models provided the motivation to introduce 
a simple compressible Sherrington-Kirkpatrick (SK) Ising spin-glass model [HI [9], and 
use standard techniques to obtain the global phase diagram in terms of pressure and 
temperature. 

The investigation of compressible Ising models has a long history flU{ [TTl [T2] . 
Mean-field and renormalization-group calculations indicate that the inclusion of elastic 
degrees of freedom may change the nature of the continuous Ising transition. An 
Ising antiferromagnet on the triangular lattice, with properly chosen elastic degrees 
of freedom, is known to provide a mechanism to explain the formation of striped 
phases p!3| ITi] . Although elastic spin models with uniform interactions have been 
much investigated, disordered compressible spin models are less explored. Nowadays 
it has been feasible to carry out detailed computer simulations for compressible models 
p!5l [T6l ITTj . including disordered compressible models, as in the work of Marshall for a 
compressible Edwards- Anderson Ising spin-glass on a two-dimensional lattice [HI [19] . 
These works have provided further motivation to study the phase diagram of the simple 
compressible SK spin-glass model. 

We treat the mean-field SK system in the pressure ensemble, which is more adequate 
to analyze the phase diagram in terms of the intensive field variables, temperature T and 
pressure p. Using the replica method, the calculation of the free energy is reduced to the 
minimization of a functional of three sets of replica variables [H]. We obtain a number 
of results in the replica-symmetric approximation. In particular, we show that the 
presence of magneto-elastic couplings introduces a simple shift in the border between 
paramagnetic and spin-glass phases, and that there appears a tricritical point along 
the para-ferromagnetic border. Also, we perform an analysis of stability of the replica- 
symmetric solution and locate the de Almeida-Thouless line [20] . In the Conclusions, we 
consider a compressible two-sublattice SK model to make contact with the experimental 
phase diagrams of mixtures of ferro and antiferroelectric hydrogen-bonded crystals. 
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2. The compressible SK model 

The compressible Sherrington-Kirkpatrick spin-glass model is defined by the 
Hamiltonian 

^ 1 

^ = - E Jij['i--l{v-vo)]SiSj-Hj2Si + ^kN(v-vof, (1) 

l<j<j<7V i=l 

where Si = ±1, with i = 1,...,N, v is the specific volume, > is a constant 
parameter, 7 is the magneto-elastic coupling, H is the external field, and A; > is 
a uniform elastic constant. As usual, {Jij} is a set of independent and identically 
distributed random variables, with suitably scaled mean values, (Jjj) = Jq/N, and 
variances, {{Jij — {Jij))'^) — J'^/N. Given a configuration {Jij}, we write the partition 
function in the pressure ensemble. 



Y 



dv exp {-(3pvN) Tr exp {-(3n) , 



(2) 



where p is the pressure, (3 = 1/ (kBT), and the trace is a sum over spin 
configurations. Performing the volume integration, and discarding irrelevant terms in 
the thermodynamic limit, we have 

^2 ■ 



p- 



fm ( pvo - — 



Y — exp 

with an effective four-spin Hamiltonian, 
'^eff = - I 1 



Trexp(-/3K//), 



jp 
k 



i<j \ i<j 



(3) 



-HY^S,. (4) 



Using a Gaussian identity, we write the partition function as 



Y — exp 
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p" 
2k 



dxexp 
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2k 
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which is in a more convenient form to be treated by replicas. 
According to the replica method. 



where 
(y") = exp 



(Iny) = lim-ln(F") , 
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with the definition 



^ k k ' 



(7) 
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Taking the average over the exchange configurations, and discarding irrelevant terms in 
the thermodynamic hmit, we have 



(F") = exp 



p 

-j3Nn ( P^^o - ^ 



+00 



exp 



2k 



xTr„exp 



2N 



a,/3 i<j 

which can be rewritten as 



i<j 



(y") = exp 



p 

-(5Nn I P^o - ^ 



I 



2k 



(9) 



Tr„exp<C^5^reM J^^r^f 
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a</3 
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(10) 



We now use Gaussian identities to introduce the new set of variables {ma}, 
associated with the magnetization, and {qa/s}, associated with the overlap between 
replicas, so that 

Y[dx'^ I Y[ dm"" I JJ dq'^^ exp [NG (x", m", g"^)] , (11) 

-00 „, J —00 „, J —00 a 



a</3 



where 

G (x", m", g"^) = -(3n (pvo - ^) - E 



2k 
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^2 , 



E (^")' 



+ In Tr exp [ J2 ^ V^'^"'^^ + /S^o E + E ) ^^2) 
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In the thermodynamic limit, the free energy comes from 

g^g{T,p, H) = lim - maxG (x", m", g'*^) , 



(13) 



where the maximum should be taken with respect to the three sets of replica variables. 



3. Replica-symmetric solution 



In the replica-symmetric solution, with Xa — x and rua — m, for all a, and qa/s — q, for 
all pairs (q!/3), we have 



P' 



■X + [l-q) z-m 



2k 
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+ / -j= exp In 2 cosh [(3 [J^q^'^z + Jo^m + H)] , (14) 

J-oo v27r V 2 / 



with 



{ = 1 + 1 + ^. (15) 



The equations of state are given by 



q = j exp (~-\ tanh^ [/3 {J^q'/'z + Mm + H)] , (16) 

m = y ^ exp ( -- J tanh [/3 ( J^g^/^^ + Um + H)] , (17) 



and 



X = 4 (1 — g j + — m . 



2 " ^ " ' 2J 
We now restrict the analysis to zero field {H = 0). 

The transition between the spin-glass and the paramagnetic phase comes from the 
expansions of equations ffTB]) and ffTS]) . 



and 



x = ^ + 0(g^). (19) 



There is a (positive) solution, g > 0, for 



l-^>0. (20) 

Introducing the notation 

t = - = ^, (21) 

and using the definition of ^, given by Eq. (fT5|) . we have the second-order boundary 
between the paramagnetic and the spin-glass phases, 

In zero field, the calculation of the para-ferromagnetic transition comes from the 
analysis of the expansion 

x = ^^ + ^m^ + 0{q'), (23) 
which can be written as 

Introducing the dimensionless and more compact notation 
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Eq. fl24l) can be rewritten as 

b_ 

Inserting this expansion into the equations of state, it is easy to show that 



e 1 



t-b 



1 + — m - - 
a Jo 



t 



(joW + 3gjo)+OK,g3/2), (26) 



from which we have the second-order para-ferromagnetic border, 



tc2 = ajo + b= 1 + 



IP 
k 



Jo + 



2k 



(27) 



The transition between the spin-glass and the ferromagnetic phases may be 
calculated from an analysis of the zero-field susceptibility. It is not difficult to show 
that this border is given by 

ta = b{l-q^)+ajoil-q), (28) 
where q comes from the equation 



+ 00 



dz 



exp 



1 tanh^ 



aq^/'^z 



t-b 1 



(29) 



In order to investigate the eventual existence of tricritical points, we write the 
expansions 



m = 

and 
q = 



Jpam 
t-b 



1 + -jom' 
a 



1 JqQ m 
^{t-bf 



2 I o --2 I o ^ ■ 4,0^ ■-! 2 

a a 



+0 (m^ q^) 



(30) 



t-l 



m 



(31) 



+6gjo~^'^^) + O (m^, q^) . 

For m = 0, the coefficients of Eq. (!3T|) indicate that there is no possibility of a tricritical 
point along the second-order border between the paramagnetic and the spin-glass phases. 
At the para-ferromagnetic border, however, it is not difficult to locate a tricritical point 

at 

3 J7J0 U! - 1) k 



Ptc 



2 (Jo' + 2) 



7 



(32) 



Let us turn to the Hessian matrix associated with the replica-symmetric solutions 
of —(3g. There is a total of thirteen distinct elements in the n — limit, seven of which 
coming from derivatives involving the variables. As in the standard SK model, the 
paramagnetic phase is always stable, and the eigenvector subspace is generated by three 
classes of vectors, which have the general form 

f m \ 

, = 1, ■ ■ ■ ,n, (33) 



{en 
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where {0"}, {e"} and {r/("^)} are column- vectors with n, n and n{n — l)/2 elements 
respectively pOj. The eigenvalue spectrum is obtained straightforwardly following the 
same steps of the rigid model. In the ordered region, the change of sign of the transverse 
eigenvalue leads to an instability. The de Almeida-Thouless line is given by 

(^^^ ={a + 2hxf exp (^-^ sech^ [{3 [a + 2hx) {Jq^'\ + Jom)] , (34) 

where H = 0, and m, q, and x, are the equations of state. 



4. Conclusions 



It is interesting to draw a typical phase diagram in terms of temperature t = ksT j J 
and the ratio jo = Jq/ J (in zero field). The dotted lines in Figure [H are the well- 
known replica-symmetric second-order boundaries, with the addition of the de Almeida- 
Thouless (AT) instability line, for the rigid SK model (with no magneto-elastic coupling). 
The location of these borders and of the AT line is shifted in the presence of a magneto- 
elastic coupling. The spin-glass and ferromagnetic regions are broadened as a result of 
a decrease of the effective energy of the system due to the magneto-elastic coupling, as 
seen in equation (jl]). The solid and dashed lines were drawn for the parameter values 
a = 1 + ''jp/k = 1.5 and b = J'~f'^/2k = 0.5. For these particular values of pressure 
and elastic parameters, the second-order boundary between the paramagnetic and the 
ferromagnetic phases, given by t = 1/2 -|- 3jo/2, ends at a tricritical point, at jo = 2, 
beyond which there is a (dashed) line of first-order transitions. These are typical results 
for the compressible SK model in the replica-symmetric approximation. 

With the exception of the well-known corrections in the spin-glass region [9], in 
analogy to the treatment of the rigid SK model, the form of the effective four-spin 
Hamiltonian in the pressure ensemble does not indicate any drastic qualitative changes 
of the phase diagrams even if we go beyond the simple replica-symmetric approximation. 

We now make contact with previous work for mixtures of ferroelectric and 
antiferroelectric hydrogen-bonded crystals of the KDP family. A disordered Ising model 
in a random field, with the inclusion of elastic degrees of freedom, has been used to 
investigate the glassy transition in Rbi_a;(NH4)2^H2As04, known as RADA [6]. Also, 
the phase diagram of Rbi_a;(NH4)a;H2P04 (RADP), in terms of temperature T and 
concentration x, has been drawn on the basis of a cluster calculation for a model of Ising 
pseudo-spins [HE]. The Rubidium crystals, RbH2P04 (RDP) and RbH2As04 (RDA), 
display a ferroelectric transition, at about 150 K, similar to the well-known ferroelectric 
phase transition in the isomorphous KH2PO4 (KDP) crystals. The Ammonium crystals, 
NH4H2ASO4 (ADA) and NH4H2PO4 (ADP), display a strong first-order antiferroelectric 
transition. At the mean-field level, the description of the T — x phase diagrams requires 
the consideration of two sublattices. We thus use a two-sublattice compressible SK 
model, given by the Hamiltonian 

^ = - 5Z J^J[l-l{v-vo)]S,S, + hN{v-Vo)\ (35) 

ieA,j€B 
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Figure 1. Comparison between the t — jo phase diagrams of the rigid (dotted hnes) 
and compressible (soUd hnes) SK models. We indicate the paramagnetic (P), spin-glass 
(SG), and ferromagnetic (F) phases. Solid lines (1) to (4) are the boundaries coming 
from the replica-symmetric solution of the compressible SK model. Solid line (5) is 
the AT instability border. There is a tricritical point separating the continuous (solid 
line) and first-order (dashed line) P-F border of the compressible SK model. 



where the sum refers to all distinct pairs of spins belonging to different sublattices, A 
and B. The random variables {Jij} are associated with a Gaussian distribution, with 



and make J'^x (1 — x), in order to mimic ferroelectric {x = 1) and antiferroelectric 

[x = 0) transitions. Two-sublattice Sherrington-Kirkpatrick models, even in the 
presence of random fields, have been studied in connection with transitions from 
antiferromagnetic to spin-glass phases [HI [221 1231 [21] . The addition of elastic degrees 
of freedom, according to the steps of the previous sections, suggests a typical phase 
diagram as sketched in Figure [21 With a convenient choice of the parameters, it is 
possible to predict an elastically induced first-order transition from the paraelectric to 
the antiferroelectric phase for x ~ 1, turning into second order for smaller values of 
concentration, and which is in qualitative agreement with the experimental findings. 
The glassy transition can still be understood within the framework of the rigid SK 
models, with the well-known instability of the replica-symmetric solution in the low- 
temperature region. It is important to remark that the inclusion of Gaussian random 
fields is not capable of explaining a first-order transition, in contrast to the predictions of 
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the proposed compressible model In fact, first-order transitions are associated 

to either a double-delta distribution of random fields [261 121] or to the consideration of 
certain discrete quadrupolar-glass and Potts-glass models [2]. 
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Figure 2. Schematic phase diagram for the two-sublattice compressiblke SK model 
for RADP in terms of temperature T and composition x. We indicate paraelectric, 
ferroelectric, antiferroelectric, and glassy regions. Solid lines represent continuous 
transitions. The dashed line is a first-order transition. 
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